Experimental growth law for bubbles in a "wet" 3D liquid foam 
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Abstract 

We used X-ray tomography to characterize the geometry of all bubbles in a liquid foam of 
average liquid fraction (f>i ~ 17% and to follow their evolution, measuring the normalized growth 
rate Q = y ^QL for 7000 bubbles. While Q does not depend only on the number of faces of a 
bubble, its average over /—faced bubbles scales as Gf ~ / — /o for large /s at all times. We discuss 
the dispersion of Q and the influence of V on Q. 

PACS numbers: 82.70.Rr, 83.80.Iz 
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Liquid foams consist of bubbles of gas separated by a continuous liquid phase occupying 
a fraction (fit of the foam's volume. Liquid foams coarsen because gas slowly diffuses through 
the liquid films from high pressure bubbles to low pressure bubbles so the high pressure 
bubbles eventually disappear. The dynamics of the coarsening differs in the two limiting 
cases of dry and infinitely wet foams. 

In a bubbly liquid, with ~ 1, bubbles are spherical and well separated. As in 
emulsions and diphasic materials, the bubbles coarsen via "Oswald ripening" as theo- 
retically described by Lifshitz, Slyozov and Wagner (LSW) [l, 2] using a mean-field approach. 

In the dry foam limit, ^<1, bubbles are polyhedral and touch each other. The gas flux 
between two neighboring bubbles is proportional to their pressure difference, hence to the 
mean curvature of the film separating them. From simple dimensional arguments 3j, the 
growth rate of an individual bubble of volume V in a 3D foam must be of the form: 

Q depends both on the shape of the bubble rescaled to unit volume and on physico-chemical 
characteristics of the liquid, gas and of the surfactant through the effective gas diffusivity 
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where H is the mean curvature of a surface element dS of the bubble. Theoretical treatment 
of coarsening in dry foams usually make an analogy with 2D foams, for which von Neumann's 
law (vNM), which is exact, states that the growth rate of the area of a bubble is proportional 
to n — 6, where n is its number of edges. This law has been validated experimentally. In 
3D, in contrast, the existence of the extra radius of curvature means that the growth rate 
of a bubble does not depend only on its number of faces WW . Thus deriving a growth law 
in 3D requires additional assumptions and results are only approximate. Attempts to relate 
the growth rates of individual bubbles or the average growth rates of classes of bubbles to 
simple geometrical properties have had a mixed success. Indeed many attempts have been 
done to find a relation between the right-hand-side (rhs) of eq. [2] and simple geometrical 
properties. Mullins [6] and several groups (see [4] and references therein) derived estimates 



for this rhs based on ideal regular bubbles and obtained Qu ~ f 1 ^ 2 for large /. Numerical 
simulations of coarsening have computed Q for large numbers of bubbles. Potts Model 
, vertex model and Surface Evolver [9] simulations all exhibit a strong correlation 
between bubbles' number of faces and growth rate with a small dispersion. We can then 
check if Mullins' model holds on average for realistic foams by evaluating the average values 
Qf =< Q >f over sets of /—faced bubbles. The data admit both Gf ~ (/ — fn) and 
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Experimental measurements of Gf using optical tomography 



15j | observed too few bubbles to distinguish (/ — / ) from (/ — / ) 



Here we investigate foam coarsening experimentally in the previously unstudied interme- 
diate wetness regime with (pi ~ 10 — 20 %, which is not known in theory nor in simulations, 
but which is widely used in most foam applications [l6]. We determine Q and Gf for several 
thousand bubbles and study the cross-over between the LSW and the vNM regimes. Recent 
progress in X-ray tomography |17j permits us to take 3D images of thousands of bubbles 
in 2 minutes, and allows us to monitor a foam for many hours. We find that, for a given 
set of /—faced bubbles, the dispersion of Q is small. Therefore we focus on Gf. While the 
shape of the distribution of the number of faces per bubble changes during the experiment 
(implying that our foam has not reached a statistically invariant scaling state), the shape of 
Gf vs. f stays constant. This quantity thus seems to be a robust parameter to characterize 
foam coarsening. Gf depends linearly on / for large f. We then study the dispersion of Q as 
a function of / and V and we finally discuss whether / or V is more relevant to determine 
the dynamics of a bubble, on the basis of the strong correlation existing between both 
quantities through Lewis' law. 



The X-ray tomography apparatus at the ID 19 beamline of the European Synchrotron 
Radiation facility (ESRF) Grenoble, France, enables us to visualize the coarsening of a wet 
liquid foam contained in a 1 cm 3 cubic volume. The spatial resolution of the images is 10 /mi 
in each direction, smaller than the smallest bubbles which appear in at least two successive 
images (i.e. any bubble smaller than 10 /mi diameter disappears in less than 2 minutes. It 
thus never enters our calculations of growth rates). Thus we can image all bubbles which we 
can analyze (by definition, we do not include in Gf bubbles which disappear within one time 
step). Obtaining adequate image contrast between liquid and gas phases is a major difficulty 
in this kind of experiment since the liquid phase occupies only a small percentage of the total 
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volume - the liquid fraction (pi [18J . In addition, the bubble walls must not move significantly 
during the time required to acquire an image otherwise topographic reconstruction of the 
foam will fail. The characteristic evolution time for gas diffusion in ~ 1 mm 3 air bubbles is 
several minutes. Our experimental acquisition time of 2 min per image enabled us to clearly 
identify and follow bubbles during their evolution. To reduce the complicating effects of 
gravitational drainage of the foam and to maintain a constant and relatively high liquid 
fraction, we used a high-viscosity liquid composed of 100 mL de-ionized water, 0.1 g sodium 
dodecyl phosphate, 0.003 g dodecanol and lg gelatine (not enough to gelify the liquid). In 
addition, we used a pump to provide a constant supply of liquid to the top of the foam 
which maintained the liquid fraction, decreasing from 20% at the top to 10% at the bottom 

0. 
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We identified bubbles in an image using a watershed method detailed in [171 ] . We measured 
the liquid fraction in a region by summing the number of pixels of that volume containing 
more than a threshold amount of liquid. <pi was thus estimated with an error Acfc i = 0.5%. 



We developed a correlation algorithm using the commercial Aphelion 19j software 



package to track bubbles from image to image, enabling us to determine their growth rates. 



Determining a bubble's number of faces in a wet foam is ambiguously defined [17|], since 
neighboring bubbles may be distant. We compared the distance separating two bubbles 
to the characteristic size of their vertices, i.e. the crossing of the bubbles' edges where 
most of the liquid accumulates. We assumed that if the distance between two bubbles was 
less than half this characteristic size, they were neighbors and shared a face. Thus small 
bubbles located at vertices may have fewer than 4 faces. We automatically discarded from 
our analysis bubbles which touched the container walls. 



48 minutes after production, the foam's liquid fraction varied from 20% at its top to 
10% at its bottom as a result of the continuous liquid injection at its top. To test the 
influence of the liquid fraction on foam dynamics, we defined four subregions corresponding 
to different heights and thus wetnesses in our column. 

for each image, we measured each bubble's number of neighbors / and growth rate Q. 
figure [1] shows Q vs. / for the 3582 bubbles in the topmost horizontal slice of the cell 
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48 min after foam creation, and the distribution of Qs for / = 12 (inset). Error bars on the 
measurement of Q result of the incertitude on the measurement of V at two successive times. 
AV was determined for each bubble by assuming that a 1 voxel thick sheet could randomly 
be added or removed from each of its faces. 
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FIG. 1: Scatterplot of Q vs. / in the topmost slice ((pi = 20%) of the experimental cell. Inset: 
histogram of Q for 12-faced bubbles. 

Q varies for different bubbles with the same fs to a much greater extent than the exper- 
imental uncertainty. Hence it is not an exact function of the unique variable /, although 
the two quantities are correlated, showing that 3D bubbles do not obey an exact equivalent 
of von Neumann's law. However for the ensemble of bubbles with a given /, the dispersion 
of Q around its average is not large, as the inset of fig. [T] shows for f=12, since one can 
clearly see the increase of Q with /: The dispersion ag coincides with the variation of Gf 
for Af « 5 (with / > 10). We therefore compute the average value Gf to investigate its 
dependence on </>/ and the structure of the foam. 

Foam structure changes during coarsening. We characterize the structure by measuring the 
distribution of bubbles' number of faces p(f). figure [2] shows p(f) at different times and for 
different values of (pi. The shape of the distribution varies significantly with time and liquid 
fraction. The average number of neighbors ranges from 10.7 ± 5.5 (</>/ = 20%) to 11.0 ± 4.5 
(4>i = 14%). Figure [2] also shows that p(f) varies significantly from 48 min to 375 min after 
foam creation, hence the structure is still far from any scale invariant regime in which p(f) 
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FIG. 2: Solid red and blue lines: p(f) for = 20%, 48 min and 375 min after foam creation. Solid 
green and black lines: p(f) for 4>i = 14% and 17%, 48 min after foam creation. 



would be constant (20]. 

Figure [3] shows Gf at two different times for different values of (pi. Despite a shift in the 
value of / at which Gf crosses 0, the general dependence of Gf vs. / is unchanged: Gf goes 
to as / —* for / < 10, and Gf is linear for / > 10. 



Since in the limit <pi — ► 1 Q should depend on V only , we plot Q vs. V in Figure HI As 
expected, both quantities correlate. 

In order to check if the number of faces describes better a bubble's behavior than its 
volume as it is expected in the dry foam limit [3] , we scatter-plot Q vs. V for three different 
sets of /—faced bubbles in fig. [51 It is clear that, at fixed /, Q increases slightly with V. 
The shift of the clouds corresponding to different /—set over the V axis can be explained 
by the strong correlation between < V > t and f, that we plot in the inset of fig. M usually 

nun 

known as Lewis' law [7|, [9|, llq . 

We experimentally observe a power law relation between < V >/ and / : Vf =< V >/~ 
f a , with a ~ 2.2 here. However this law does not account for the slight increase of Q vs. V 
for / = 16. 
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FIG. 3: Blue and red points Gf vs. / at times 48min and 375 min after the making of the foam 
for (pi = 20%. The standard deviation of Q around Gf for each / is pictured by the bars. Green 
and black points : Gf vs. / at time 48 min for = 14% and <pi = 17%. 

We compared the dispersion of the data points for scatter plots of Q vs. / (fig. d]) 
and vs. V (fig. HJ). for / = 12 (476 bubbles) we obtained a standard deviation for Q of 
AQ = 10.5 /im 2 .s _1 . On the basis of the empirical Lewis' law Vf ~ f 2 ' 2 we plotted on fig. [5l 
we compute a AV^ corresponding to Af = 1 and find AQ = 12.1 /im 2 .s _1 around Vi2- This 
small difference does not allow to conclude that one parameter is better than the other. 
Future experiments using foams of very different structures (bidisperse, monodisperse, 
scaling state...) and liquid fractions may allow to obtain various relations between / and 
V, and thus to determine which parameter, V or /, influences Q most significantly. 



Our experiments studied a moderately wet foam (14% < 0; < 20%), for the dynamics of 
which we have no theoretical predictions. Q does not depend uniquely on / and Gf disagrees 
with mean-field theories for both wet and dry foams. Gf increases linearly with / for / > 10 
as Glazier numerically observed using Potts' simulations . Thomas et al [7] have confirmed 
Glazier's results and extended them to additional initial bubble size distributions. As in dry 
foam simulations, Gf increases to zero as / decreases for small fs js, 3]. Let us notice that 
there are still some spurious bubbles due to the noise in image segmentation. Some of them 
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FIG. 4: Growth rate vs. V, volume of the bubbles inside the topmost slice (4>i = 20 %) in the 
48 mm old foam. 
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FIG. 5: Scatterplot of the growth rate Q vs. V for 8 (red), 12 (blue) and 16-faced (green) bubbles. 
Inset: < V >/ as a function of / for (f>i = 14, 17 and 20%. 
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are artificially correlated, but their number has been evaluated to be very small, and their 
influence on the results is poor as was checked by introducing a filter on the size of the 
bubbles taken into account. 

In contrast to 2D, the dispersion of Q around Gf is intrinsic to the foam, as / does not 
suffice to describe the dynamics of a bubble. In addition we have checked the stability of 
the linearity of Gf for early stages of the evolution of the foam. The foams contains too few 
bubbles at late stages to compute statistically relevant GfS. It will be of great interest to 
observe Gf in later stages, in particular during a possible scaling regime. 
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